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Q ' Angular momentum transport by internal magnetic fields is an important in- 

^ ■ gredient for stellar interior models. In this paper we critically examine the basic 

heuristic assumptions in the model of the Tayler-Spruit dynamo, which describes 

^ how a pinch-type instability of a toroidal magnetic field in differentially rotating 

>■ ■ stellar radiative zones may result in large-scale fluid motion. We agree with prior 



published work both on the existence of the instability and its nearly horizontal 



in 

O ' geometry for perturbations. However, the approximations in the original Ache- 

Q ■ son dispersion relation are valid only for small length scales, and we disagree that 

^ ■ the dispersion relation can be extrapolated to horizontal length scales of order 

the radius of the star. We contend that dynamical effects, in particular angular 

P^' momentum conservation, limit the maximum horizontal length scale. We there- 

O ' fore present transport coefficients for chemical mixing and angular momentum 

5— ( ' 

-t— » ' redistribution by magnetic torques that are significantly different from previous 

ci ■ published values. The new magnetic viscosity is reduced by 2 to 3 orders of mag- 

nitude compared to the old one, and we find that magnetic angular momentum 
^ ' transport by this mechanism is very sensitive to gradients in the mean molecular 

j^ ■ weight. The revised coefficients are more compatible with empirical constraints 

on the timescale of core-envelope coupling in young stars than the previous ones. 
However, solar models including only this mechanism possess a rapidly rotat- 
ing core, in contradiction with helioseismic data. Previous studies had found 
strong core-envelope coupling, both for solar models and for the cores of massive 
evolved stars. We conclude that the Tayler-Spruit mechanism may be important 
for envelope angular momentum transport, but that some other process must be 
responsible for efficient spindown of stellar cores. 
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1. Introduction 

Mixing driven by stellar rotation is important for a variety of problems in stellar struc- 
ture and evolution. Hydrodynamical mechanisms, such as meridional circulation and shear 
instabilities, induce mixing that is consistent with observational data (see Pinsonneault 1997; 
Maeder & Meynet 2000 for reviews on rotational mixing in low and high mass stars, respec- 
tively). However, the evolution of a rotating star depends on its internal angular momentum 
distribution, and existing theoretical models have difficulty matching empirical constraints 
on the timescale for transport in stellar radiative interiors. The slow rotation of the solar 
core (Tomczyk et al. 1995) and newly born pulsars (e.g., Ott et al. 2005) indicates that 
angular momentum is transferred more effectively than would be expected from hydrody- 
namical mechanisms alone. To complicate matters further, solid body rotation is also not 
an appropriate approximation for important phases of stellar evolution. Both the spindown 
of young cluster stars (Stauffer & Hartmann 1987; Keppens et al. 1995; Bouvier et al. 1997; 
Krishnamurthi et al. 1997) and the survival of rapid rotation in old horizontal branch stars 
(Peterson 1983; Behr et al. 2003; Sills & Pinsonneault 2000) require that internal differential 
rotation persist for timescales of order 20 Myr to 100 Myr. 

There are two other major classes of mechanisms that could be important agents for 
angular momentum transport: magnetic fields and waves (g-modes). Both waves and mag- 
netic fields have the property that they generally transmit angular momentum much more 
effectively than they induce mixing. As a result, the efficiency of these processes is most 
easily tested by their impact on angular momentum evolution. In this paper we focus on 
magnetic angular momentum transport via the Tayler-Spruit mechanism. 

Elaborating upon results of earlier investigations of instabilities of toroidal magnetic 
fields in stars conducted by Tayler and others (e.g., Tayler 1973; Markey & Tayler 1973, 
1974; Acheson 1978; Pitts & Tayler 1985), Spruit (1999) has proposed a new magnetohydro- 
dynamical mode of angular momentum transport in radiative zones of differentially rotating 
stars. Fluid elements experience large-scale horizontal displacements caused by an unstable 
configuration of the toroidal magnetic field (one consisting of stacks of rings concentric with 
the rotation axis). Small-scale vertical displacements of fiuid elements are coupled to the 
horizontal motions, which can cause both mild mixing and much more effective momentum 
transport. 

Spruit's key idea is that no initial toroidal magnetic field is actually required to drive 
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the instability and mixing because the unstable field configuration can be generated and 
maintained by differential rotation itself in a process similar to convective dynamo. The 
dynamo cycle consists of two consecutive steps: first, a poloidal field Br is generated by 
the vertical displacements of the unstable toroidal field; second, the new poloidal field is 
stretched into a toroidal field B^ by differential rotation (for more details, see §2). 

Following Spruit (1999, 2002), who considered two opposite limiting cases in which 
the thermal diffusivity K is either negligible or dominating, Maeder & Meynet (2004) have 
derived equations for the transport by the Tayler-Spruit dynamo applicable to a general case. 
These equations as well as Spruit's original equations have already been applied to study the 
angular momentum transport by magnetic torques, which are proportional to the product 
B^B^, in massive stars on the main sequence (MS) (Maeder & Meynet 2003, 2004, 2005) and 
during their entire evolution toward the onset of iron-core collapse (Heger et al. 2005; Woosley 
& Heger 2005). While the Geneva group has shown that the Tayler-Spruit dynamo can easily 
enforce the solid-body rotation of massive MS stars, Heger and co-workers have established 
that appropriate mass-loss rates coupled with the magnetic torques generated by the Tayler- 
Spruit dynamo can spin down pulsars to the observed rotation rates. However, the same 
angular momentum transport model has failed to explain gamma-ray bursts (Woosley & 
Heger 2004). Finally, Eggenberger et al. (2005) have recently demonstrated that the Tayler- 
Spruit dynamo leads to a nearly flat solar rotation profile consistent with that revealed by 
the helioseismic data (e.g. Couvidat et al. 2003). The important question of stellar spindown 
was not directly addressed in this theoretical work. 

In this paper, we derive new transport coefficients for the Tayler-Spruit dynamo that 
have quite different magnitudes and dependences on the /i-gradient (§2 and §4). We review 
prior work in §2. The diffusion coefficients are obtained as products of the growth rate 
and the characteristic vertical length scale of the instability. Spruit (1999) used a local 
dispersion relation to infer the ratio of vertical to horizontal length scales. To convert 
this ratio to an absolute vertical length scale he assigned the maximum possible horizontal 
length scale r. We contend that the assumptions in the dispersion relation break down 
for large length scales, and that the maximum effective horizontal length scale cannot be 
obtained by the method employed by Spruit. In §3, we give some geometric, energy and 
dynamical arguments that support reduced length scales and we present revised diffusivity 
estimates. The original prescription predicts extremely rapid core-envelope coupling for 
low mass stars (in contradiction with observations), but also a flat solar rotation curve 
(in agreement with observations). The situation is reversed with our diffusion coefficients; 
the spindown timescale is closer to that inferred from young open cluster stars, but a solar 
model computed with the new equations possesses strong differential rotation in its radiative 
core, in conflict with helioseismology. In § 4, we discuss the observational consequence. We 
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therefore conclude that some other magnetic mechanism, or waves, must be responsible for 
angular momentum coupling between deep stellar cores and envelopes. 



2. Revised Heuristic Assumptions and Transport Coefficients 

The diffusion coefficient estimates for the Tayler- Spruit mechanism were derived on a 
heuristic level, and the results were confirmed by an appeal to a restricted form of the full 
Acheson (1978) dispersion relation. We begin this section by summarizing the overall argu- 
ments in Spruit (1999, 2002), noting the places where we differ with him. In particular, we 
claim that the results from the dispersion relation can only be applied to small horizontal 
length scales, therefore they cannot be used to justify a characteristic length scale for hori- 
zontal motions of order the radius of the star. In addition, the presence of curvature terms 
and variations in physical quantities such as the rotation rate during the motions will make 
both our estimates and those of Spruit upper limits. 

The basic parameters of the Tayler-Spruit dynamo include the instability's growth rate 
a;, the magnetic diffusivity r], buoyancy (or Brunt- Vaisala) frequency N, angular velocity Q 
and its gradient q = [dlnfl/dlnr) in a stellar radiative zone. The vertical displacement Ir 
of a disturbed fluid element has an upper hmit 



(appropriate expressions for ui and iV^ are provided below). The displacement is constrained 
by the requirement that its kinetic energy must exceed the work that has to be done against 
the Archimedes force. At the heuristic level. Spruit justified this equation in terms of energy 
conservation. He balanced the magnetic energy gain (from motions perpendicular to the 
rotation axis) with the work done against the buoyancy force (from motions along the radial 
direction). This is an upper limit, as it assumes that all of the magnetic energy can be 
converted into vertical motions. He justified this in the dispersion relation by adopting a 
simplified version of the full Acheson (1978) equation valid only near the rotation axis. In 
this case, the horizontal and vertical directions correspond to motions perpendicular and 
along the rotation axis respectively. The ratio is maximized at the pole, which is another 
way in which this assumption serves as an upper limit rather than a bound. One can obtain 
a constraint (equation A21 in Spruit 1999) on the ratio of the horizontal and vertical length 
scales. In cylindrical polar coordinates {w,ip,z), used by Acheson and Spruit, this has the 
general form l^jUJ^ > I'^N'^ for an m = 1 mode and Spruit's a = 1. In the dispersion relation 
that was used, there is no bound on the energy gain from horizontal motions; Spruit set 



-5- 



the horizontal length scale to be the radius of the star, which is certainly a maximum. 
We will argue in § 3 that the large scale horizontal fluid motions initiated by the Tayler 
instability are restricted by the Coriolis force within bounds of order /^ <^ r, and that the 
basic assumptions in the dispersion equation, for instance those that both /^ and l^ should 
be much smaller than r (see § 2 in Acheson 1978), imply that this equation cannot be used 
to extract information about the absolute maximum horizontal displacement. In particular, 
the subtle assumption that motions perpendicular to the rotation axis do no work against 
the buoyancy force is valid only for small perturbations close to the rotation axis. 

There is also a lower limit 



I'r > -, (2) 

UJ 

where r] is the microscopic magnetic diffusivity (i] = 7 x 10^^ In AT^^/^ cm^ s^^, where InA ^ 
5-10 is the Coulomb logarithm). This bound is determined by the length scale at which 
the magnetic diffusion dissolves any new poloidal field faster than it is produced by the 
instability. According to Pitts & Tayler (1985), the instability's growth rate takes on a value 



u = uja = J-. — ^75—7 when il <^ lja, (3) 



(47rp)i/2r 



and 



^2 
cu = -y^, when n > cua, (4) 

where u^ is the Alfven frequency, p and r are the local density and radius. 

We are mainly interested in the case of fast rotation, when Q ^ uja- This condition is 
satisfied even for the relatively slow solar rotation. In this case, the Coriolis force reduces 
the instability's growth rate considerably by a factor of uja/^ ^ 1 (Pitts & Tayler 1985). 

In § 3, we will argue that one should employ the same value of u, either (3) or (4), when 
estimating the heuristic limits (1) and (2). Below, we present revised diffusivities obtained 
under this basis. However, Spruit (1999) and, after him, other researchers (e.g., Maeder & 
Meynet 2004) have substituted the plain Alfven frequency (3) into inequality (1) but the 
reduced growth rate (4) into inequality (2) and all other equations. Having done that, the 
appropriate heuristic expressions for the effective magnetic diffusivity rje ^ l^oj and viscosity 
fe = BrB^/{ATipqVL) are obtained quite straightforwardly. We therefore repeat the steps in 
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their derivation and note the places where our different assumptions yield different conclu- 
sions. When deriving them, we will adhere to the following three important assumptions 
made by Spruit (1999) and Maeder & Meynet (2004). First, they have considered the case 
of marginal stability, i.e. /^ has been set equal to its upper limit (1). 

From the magnetic induction equations for Br and -B^ (see below) it follows that the 
joint operation of the Tayler-Spruit instability and differential rotation can reach a stationary 
regime, in which the poloidal and toroidal fields maintain each other, only if 



Indeed, let us denote B = i^^e,. a weak poloidal magnetic field (here, we use the spherical 
polar coordinates). If, for the sake of simplicity, we only consider gas motion in the equatorial 
plane then its velocity due to rotation is v = flre^p. Let us also neglect the magnetic 
diffusivity rj for a while. Under these assumptions, the magnetic induction equation 

dB 



:_ = V X (v X B) + r]V^B 



is reduced to 



or, after choosing a frame of reference rotating with the angular velocity Q, 

^ = r|^B. = SJ,B„. (7) 

where q = (dlnfl/dlnr). This equation shows how the differential rotation winds up a 
toroidal field B^e^ by stretching the poloidal field lines around the rotation axis. However, 
the azimuthal field is subject to the Tayler-Spruit instability at any field strength (Pitts & 
Tayler 1985; Spruit 1999). This non-axisymmetric pinch-type instability, whose driving force 
is the magnetic pressure B^/Sn, causes the concentric field lines to slip sideways a horizontal 
distance /h accompanied by a radial displacement Ir- From the Acheson equation. Spruit 
derived the estimate Ir/lh ^ ^a/N. 

Like the differential rotation generates the azimuthal magnetic field when operating 
on the seed poloidal field (eq. 6), the radial displacement /^ with a characteristic velocity 
f r ~ /r c^ produces a weak poloidal magnetic field BrGj. at the expense of the existing toroidal 
field: 
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After that, the differential rotation can operate on the freshly generated poloidal field to 
wind it up into a new azimuthal field, and so on. Hence, a coordinated action between the 
differential rotation and the radial displacements caused by the Tayler-Spruit instability can 
work as a dynamo, provided that some conditions leading to a stationary regime are realized. 
The stationary regime is achieved when the magnitude of the poloidal field Br generated by 
the instability during its growth time t ^^ u"^, as described by equation (8), coincides with 
the magnitude of the field Br in equation (7) whose lines can be stretched by the differential 
rotation, during the same time r, into an azimuthal field of the same magnitude B^ as the 
one that started this dynamo loop. Quantitatively, these requirements are expressed as 

Br^r^^ ^fi^ ^ ^r^ = ^-^TQqBr, (9) 

ot r r ot r 

which results in equation (5). 

Third, one must take into account the reduction of the stable thermal stratification of 
the radiative zone caused by heat losses from surfaces of the disturbed fluid elements during 
their vertical displacements. Following Maeder & Meynet (2004), for rj^ 'ti K the reduced 
buoyancy frequency can be expressed as 






N= ^;i:^ N^ + Nl (10) 

where 



iV| = |-(Vad-Vrad), and Nl=g^ 
Hp 



d\nji 



dr 

are the T- and /i-component of the square of the buoyancy frequency in the absence of mixing. 
Here, Vrad and Vad are the radiative and adiabatic temperature gradients (logarithmic and 
with respect to pressure), Hp is the pressure scale-height, g is the local gravity, ji is the 
mean molecular weight, and 

is the thermal diffusivity with k and Cp representing the opacity and the specific heat at con- 
stant pressure, respectively. The quantities 5 = — ((91np/91nT)p and ip = (9 In p/9 In /i)p^ 
are determined by the equation of state. The value of the constant C in equation (10) de- 
pends in general on the assumed geometry. Following Maeder & Meynet (2004), we adopt 
C = 2. 

The effective magnetic viscosity can be expressed in terms of r]c' 

-^) » r/e. (12) 



The last formula has been derived using magnetic induction equations for Br and B^p (see 
eqs. 7 and 8) in a way similar to that described by Spruit (1999). 

To summarize, in general we have 9 independent equations 






Br = iB,, iV^ = r^^iVl + iV2 ,^, = /2^„ 



AirpqQ ' 
for 9 unknown quantities 

/r, Ua, Br, B^, N, r]e, l^e, (^1, and U2. 

For the case of fast rotation, Spruit has taken ui = ua and uj2 = uj\/Q. In Appendix 
A to their paper, Maeder & Meynet (2005) have assumed that u;i = C(J2 = ^a? which is 
correct for the case of slow rotation considered there. Unlike Spruit (1999), we have chosen 
uji = LJ2 = u\/Q. In § 3, we use arguments from energetics, dynamics, and a discussion of 
the limitations of the usage of the dispersion relation to justify this choice. 

Assigning a value of the upper limit (1) to /^ and comparing it with the expression for /^ 
from (5), we obtain A^ = qQ. After substituting this into (10), we find the effective magnetic 
diffusivity 

In order to compare our revised transport coefficients for the Tayler-Spruit dynamo 
with those used by Maeder & Meynet (2004), the latter are summarized below in a slightly 
modified form: 

r]e = aKy^, and u^ = a K ——Y^y"^ . (14) 

Here, y is a solution of the 4th order algebraic equation 

ay^-ay^ + (3y-2 = 0, (15) 

where 

fi^ q'^ N'^ 



-9- 



are dimensionless coefficients. 

Originally, Spruit (1999, 2002) only considered the two limiting cases: N^ ^ A*"^ and 
Nt <^ A^;^. In the ffist case, we can neglect the /3-term in equation (15) and re-write it as 
y^iy — 1) ^ 2/q;, where 

i?0 ) V 10-5 j \K )\ N^ 



„.«4xio-M^ #. m^).^ (17) 



The normalizations in the last equation assume that all quantities are expressed in cgs units. 
In the Sun's radiative core (at least, in layers located at some distance from both its center 
and the bottom of its convective envelope), r < Rq, Vt < 10"^, 10^ < Ji' < 10'', Nt > 10"^, 
and g <^ 1. Hence, there we have a <^ 1, which results in y ^ (2/q;)^/^. Using this 
approximate solution, we estimate 

r/e ^ 23/V/^ir ^ 7.89 X lOM -^ ( -^ ( — \ (^ g, cm^s-^ (18) 



Rq) Vio-v Viov V Nr 



Except for the numerical coefficient and parameter normalizations, this equation appears to 
be identical with eq. (43) from Spruit (2002) and eq. (23) from Maeder & Meynet (2004). 

For the opposite case of A^^ ^ Nt-, we have /5 ~ 2. Therefore, eq. (15) can be re-written 
as {ay'^ + 2)(y — 1) ^ 0. Since a, with Nt replaced by A^^, still remains a small quantity, 
y ^ 1 is now a solution. In this case, 

r/e ^ ai^ ^ 4.84 X 10^ (^-^y (^^ {^-^^\\ cm' s'\ (19) 

This equation coincides with eq. (42) from Spruit (2002) and eq. (21) from Maeder & 
Meynet (2004). 



„. « 2.00 X 10= f 4) (-^X ('-^)\'. '^n,^-'. (20) 



In a region where A"^ <^ Nt, our magnetic diffusivity (13) is estimated as 

Its values are smaller than those given by the old prescription (18) by a factor of ~ 10^. In 
the same region, the magnetic field strengths (in Gauss) are 

B,«9.88xl0V"^(^j (^j (^j (^j f", (21) 
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in all three cases considered above (eqs. 18-20), the effective magnetic viscosity and 
diffusivity are related to each other by equation (12), which gives 

/ \ ^/^ / O \ V3 / \ 2/3 

z/e = 3.65 X loM -^ -^ P »r/e. (23) 



RqJ Vio / V? 

In our case, as has been shown, N = Qq. Interestingly, the same constraint on the 
buoyancy frequency has also been obtained by Maeder & Meynet (2005) in Appendix to 
their paper, where they considered the alternative case of slow rotation, i7 ^ cja- This 
coincidence is not surprising because they have substituted the plain Alfven frequency into 
the lower limit (2) as well, as prescribed by equation (3) for slow rotation. This also explains 
why our revised effective diffusivity (13) coincides with theirs. However, our expression for 
the magnetic viscosity (12), which controls the transport of angular momentum, is different 
from that derived by Maeder & Meynet (2005) (they got u^ = f]e) for slow rotation because 
we consider the case of fast rotation with uj = uj\/VL in all equations. If we had examined 
the slow rotation case, we would recover Maeder & Meynet's result in that limit. Finally, it 
can easily be shown that our revised equations pass the consistency check demanding that 
the rate of magnetic energy production is equal to the rate of the dissipation of rotational 
energy by the viscosity z/g. 



3. Arguments Supporting Our Heuristic Assumptions 

We begin the justification of our heuristic assumptions by critically analyzing the validity 
of the form of the dispersion relation employed by Spruit (1999)^. In this initial section we 
demonstrate that the approximations in the dispersion relation at the rotation axis are valid 
only for small horizontal length scales and cannot be extended to a length scale comparable 
to that of the star, as was assumed by Spruit. The dispersion relation sets only the ratio 
of length scales, and can only rigorously be used for perturbations. Our strongest result is 
therefore that the coefficients presented by Spruit must be treated as upper bounds. 

We then follow up with a discussion of the dynamics to obtain other information about 
behavior for larger displacements. We argue that rotation limits the maximum horizontal 



^In the original version of our paper (astro- ph/0604045), our mathematical justification based on a re- 
analysis of the Acheson dispersion relation contained subtle flaws. Our conditions (24) were actually sufficient 
but not necessary conditions for the development of the Tayler-Spruit instability (we are grateful to Tony 
Piro for bringing our attention to this). Although now we agree with the estimate (A21) obtained by Spruit 
(1999), we disagree with his implementation of it, particularly with Spruit's choice of the horizontal length 
scale. 
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length scale with the same damping factor as present for the small scales where diffusion 
becomes important. The energy constraints from the magnetic induction equations are then 
considered, and we conclude with some general comments. 



3.1. Geometric Arguments 

In the form of the Acheson dispersion equation used by Spruit (1999), /h actually stands 
for a distance from the rotation axis, i.e. /h = Iw in the cylindrical polar coordinates {w, ip, z), 
while Ir = /z is a displacement along the rotation axis. Let R denote the radius of a magnetic 
ring that initially is concentric with the rotation axis and lies on a sphere of the radius r. 
The Tayler-Spruit instability will cause a leading edge of the ring to move a distance /^ 
perpendicular and a distance Iz parallel to the rotation axis (Fig. 1). As a result, a radial 
displacement of the leading edge will be 



(24) 



Because Spruit's form of the dispersion relation is defined at the poles, motions in the 
w direction are assumed to be perpendicular to gravity; this is not in general true for a 
spherical geometry. We can check on the maximum length scale where this approximation 
breaks down as follows. Assume that we start at the rotation axis {R = 0) and have /^ = 0. 
We consider that the assumptions in the dispersion relation break down when the energy 
loss from motions in the w direction {I'^N'^) equals the energy gained by the magnetic field 
{r'^uj\). Under these assumptions, equation (24) can be re- written as 




'' 'At) -' = i«^- 



from which it immediately follows that l^ < ^/3{uja/N) r <^ r. Indeed, in the deep solar 
core N ^ N^ > Nt ~ 10~^, while a;A ^ ^ ~ 10~®, hence /^ < 0.05 r <^ r. For longer 
horizontal length scales than this, the neglect of curvature terms in the dispersion relation 
cannot be justified; motion that is assumed to gain energy in the dispersion relation actually 
loses energy in a spherical star. This is consistent with Acheson's discussion of the dispersion 
relation, where he notes that it should be used only locally and not applied globally. We 
therefore conclude that the assumption /«, ~ /h ~ ^ cannot be obtained rigorously from the 
dispersion relation, and that other considerations (such as energy and momentum conser- 
vation) should be used instead to determine the maximum possible length scales. However, 
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the relationship between the horizontal and vertical length scales derived from the dispersion 
relation is valid locally near the pole. 

Another way of looking at this problem is to consider the consequences of Spruit's 
assumption that /^ ~ /h ~ ''" in equation (24). If we take into account that R cannot be 
less than ly^ (otherwise, the ring will intersect the rotation axis which is not allowed in the 
Tayler-Spruit instability) and that Iz < Iw{uja/N) <^ ly, (this is a correct form of eq. 25) 
then we will find that /j. ~ r as well. This is in contradiction with Spruit's estimate of 
Ir < r{ujx/N) <^ r. Equation (24) shows that Ir can be much less than r only if /^ is much 
less than r as well. Thus, this simple analysis demonstrates that: (i) the assumption of 
//i ~ r is ruled out by the geometry of the problem, and (ii) even our revised diffusivities 
obtained under the assumption that W^ r {oj^/^) ^ r (see below) may give overestimated 
values because, like Spruit, we did not discriminate between /h and Z^, and between 1^ and Iz- 



3.2. The Coriolis Force 

Our revised magnetic diffusion coefficient and viscosity differ from those obtained by 
Spruit (1999) only because we have used a different (reduced by the factor of oj^/VL) upper 
limits for the horizontal and radial wavelengths of the unstable fluid displacement in (1). In 
this section, we show that our choice of l^ is physically more consistent than that made by 
Spruit as we take into account a reduction of the horizontal wavelength /^ by the Coriolis 
force. 

To prove that his heuristic derivations yielded correct upper and lower limits for Z^, 
Spruit analysed the Acheson (1978) dispersion relation using a few simplifications. The 
original Acheson equation includes all relevant physical processes in an electrically conducting 
spherical body rotating in the presence of a toroidal magnetic field and its own gravitational 
field. The analysis of the dispersion relation in the limit of fi ^ cja has led Spruit to the 
conclusion that in the case of -ft' = 0, which is applicable when A^^ ^ Nt, sufficient and 
sufficiently accurate necessary conditions for oscillatory modes of the Tayler-Spruit instability 
are both met if 

and 

Whereas the second condition is already identical with the heuristic inequality (2), the first 
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condition would coincide with the inequahty (1) only if the horizontal wavelength lii = r 
(assuming, of course, that c<j = cja in (1), as Spruit did). 

The important point we want to emphasize here is that the condition (25), derived from 
the Acheson equation, does not give a direct estimate of the radial length scale Ij.. Instead, 
it says that the maximum possible value of /,. is proportional to l^ with a small coefficient 
uja/N. In order to get a final estimate for Ir, we have to make a physically motivated choice 
of /h- Spruit chose "the longest possible horizontal wavelength" Z^ ~ f^- But his choice 
is based on a purely geometric constraint which can be applied only if there are no other 
constraints that bound the unstable displacement along the horizontal surface. We argue 
that Spruit's choice is only valid in the absence of rotation (see, however, §3.1). 

In a rotating star (here, we assume the uniform rotation, like Spruit did when he analysed 
the Acheson equation), the Coriolis force 2v x J7 can considerably reduce the horizontal 
length scale of fiuid motions. To demonstrate this, we assume that the leading edge of 
the ring has the initial velocity v = {va, 0} directed off the rotation axis along the X axis 
(Fig. 3). Here, va = tua is the Alfven velocity, and the coordinate plane XOY is chosen to 
be perpendicular to the rotation axis. Our assumption is consistent with the development 
of the Tayler instability. Having solved the momentum equation, we find that the motion of 
the ring under the action of the Coriolis force is described by the following equations: 

X = Rcos9 + r — sin(2Qt), 

y = Rsme + r^[cos{2nt)-l], (27) 

where R is the radius of the ring, and the angle < ^ < 27r specifies a position of a point 
on the ring. In the limit of fi ^ (no rotation), these equations are transformed into 
X -^ Rcos6 + rujAt and y -^ RsinO, i.e. they describe a displacement of the ring with 
the constant velocity va along the X axis. However, when Q ^ the Coriolis force will 
defiect the ring from this rectilinear motion. For example, for a particular choice of r = 10, 
i? = 2 (in arbitrary units), and the ratio coa/^ = 0.1, consecutive positions (for the phases 
2Qt = 0, tt/2, n, and 37r/2) of the ring are plotted in Fig. 3. From equations (27) it follows 
that, in the case of uja <^ fi, the Coriolis force restricts the ring motion within a distance 
of /h ~ ''"(^a/^) around its original location. In the limit of uja/^ -^ 0, equations (27) are 
reduced to a; ^ Rcos6 and y — > i?sin^, i.e. in this extreme case /h -^ 0, and the ring cannot 
move horizontally at all. 

It is interesting that a similar restriction imposed by the Coriolis force has been studied 
by Choudhuri & Oilman (1987) in relation to the axisymmetric rise of a magnetic ring in 
the solar convective zone caused by the buoyancy instability. Like in our case, although the 
Acheson equation gives the same reduced growth rate u ~ uj\/Q for the buoyancy instability. 



-14- 



it does not predict how the magnetic ring is deflected by the Coriohs force on the large length 
scales. To find this out, one needs to consider the dynamics of large scale displacements. 

What happens if Q is not held constant during the non-axisymmetric oscillations? For 
a ring of size R perturbed by a distance /^ from the rotation axis, angular momentum 
conservation implies 

d2 

^ = ^ . 

'R' + il 

There will be an oscillating angular velocity gradient produced by this with a characteristic 
vertical length scale /^ of magnitude 



dz h hR' + lV 

hence 

z dQ zIIj 

Since Iw/h ^ 1, this implies a large q unless zlyj/{E^ + O ^ I5 which corresponds to small 
motions of large rings. A back reaction would therefore be produced, damping the motions, 
unless VL ~ constant were maintained, hence this is a reasonable assumption. 

Furthermore, taking into account the fact that the Coriolis force induces the azimuthal 
motion, the assumption of the star's uniform rotation made by Spruit seems inconsistent. 
The natural question arises — what causes the corotation of the fiuid element with its 
surroundings when it moves off the rotation axis along the horizontal surface or, in other 
words, what force compensates the Coriolis force? In Spruit's derivation, this cannot be a 
viscous friction because viscosity is neglected in the version of the Acheson equation used by 
him. A relevant example from the Earth's atmosphere science is the geostrophic wind. This 
wind can flow along a meridian only thanks to the presence of a horizontal pressure gradient 
that balances the side action of the Coriolis force. In stars, there is no such gradient on the 
equipotential surfaces. Rot at ion- driven shear instabilities in stars may result in a strongly 
anisotropic turbulent viscosity with its horizontal component strongly dominating over the 
vertical component (Zahn 1992). So, in fact, it is rather the horizontal viscous friction, not 
taken into account by Spruit, that may convert the azimuthal kinetic energy of the fluid 
element into the kinetic energy of the horizontal turbulence. 



3.3. Support from the Magnetic Induction Equations 

Our result is derived by limiting the horizontal length scale in the presence of rotation. 
The Coriolis force bounds /h (or, to be more precise, /^) by the limit of order r{ujp^/Q), 
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which ensures the angular momentum conservation. In this section, we demonstrate that 
our reduced upper hmit on /h, combined with the hmit (25) on the ratio of the radial 
and horizontal length scales obtained by Spruit (1999), is consistent with the requirement 
that the excess kinetic energy in differential rotation — the only source of energy in the 
problem — always exceeds the work needed to be done against the buoyancy during the 
radial displacement. 

Assuming that all of the requirements (9) are satisfied, we can write 

Bl ^ Q'q'By, (28) 

where r ~ u^^ = VL/uj\. On the other hand, 

Bl^^f^Bl = A7,pllul. (29) 

In the last equality, we used the definition of cja (eq. 3). Combining equations (28-29), we 
obtain 

g4p.^.i.lpnv<j(^)'\ (30) 

Earlier we found that A^ = Vtq. There is a natural energetic argument for this (Maeder 
2006, private communication). The product \pVl?'q^ll measures an excess of kinetic energy 
(per unit volume) in the differential rotation between spherical shells separated by a radial 
distance Ir. This is the only source of energy available to do the work ^pl^N^ against the 
Archimedes force to produce mixing on the length scale Ir in the radial direction. Without 
this mixing, no radial magnetic field will be generated by the magnetic induction. This 
conclusion is especially evident when we consider a deep radiative core of the star where 
N"^ ^ A^. Here, mixing itself cannot affect the value of A^^. If fi^g^ is less than A^ then 
the differential rotation simply does not possess a sufficient amount of the excess kinetic 
energy to overcome the stable thermal stratification. It does not matter what sophisticated 
mechanisms we may invent to transform this energy into other forms. In the end, it will 
still have to exceed the work |p/^A^ and, by the energy conservation law, it will still be 
equal to ^pVL^q^t^ at most. This is why our revised magnetic diffusivities vanish whenever 
Vt^q^ < A"^. Substituting the inequality Vl?q^ > N"^ into (30), we immediately get our 
estimate of Ir < ruj\/NVL. 

It should be noted that eq. (30) is in concordance with our conclusion that the Coriolis 
force reduces the horizontal length scale by the factor of o'a/^, provided that ^pfi'^q'^l'^ = 



^pllN"^. Indeed, under this assumption eq. (30) can be written as 



2 



l_o2„2i2 -'-„r2Ar2 -*- „ »,2, .2 /'^A 



l^P^'<iX = l^PnN' = l^pr'<[^) ■ (31) 
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After the subsitution of /^ = hi{^A/N) into the mid term and the comparison of it with 
the last term, we will immediately find that /h = r{ujA/^)- This means that the latter 
relation is actually equivalent (in the mathematical sense) to the energy balance equation 
^pVt^q^ll = ^pl'^N'^. Equation (31) shows that only a small fraction of the energy of the 
toroidal magnetic field is used to do the work against the buoyancy. On the other hand, 
if we had followed Spruit's assumption that jpr'^uj^ = jpl'^N^ then we would have found 
from eq. (30) that \pVl?q^ll <^ ^pl^N"^. In that case, some extra energy source, besides the 
differential rotation, would have been needed to drive the Tayler-Spruit dynamo. 



3.4. A Common Sense Argument 

Equations (18) and (23) show that, in the limit of N^ ^ A*"^, Spruit's magnetic viscosity 
does not depend on the shear q at all. This is difficult to understand conceptually because 
it is the differential rotation that makes the whole Tayler-Spruit mechanism work. In our 
revised prescription, the dependence of diffusivities on q holds in all regimes. 



4. Comparison with Observations 

The original Spruit (2002) diffusion coefficients for angular momentum transport were 
both large and independent of the angular velocity gradient in the absence of p gradients, 
and mildly sensitive to p gradients. The Maeder & Meynet (2004) coefficients have been 
derived in a more self-consistent fashion, but exhibit the same global properties. Our revised 
prescription has fundamentally different dependencies on both Vt and p gradients, which 
has an important impact on the behavior of stellar models. One important test of the 
impact of the Tayler-Spruit mechanism is the inferred solar rotation curve, which we estimate 
below. The older diffusion coefficient estimates were consistent with the solar core rotation 
(Eggenberger et al. 2005), and the revised coefficients predict core rotation too high to be 
compatible with seismology. 

However, there are two other important observational constraints on angular momentum 
transport in low mass stars that are also important: the spindown of young cluster stars and 
the preservation of sufficient angular momentum in red giant cores to explain the rapid 
rotation of horizontal branch stars. In both cases the observational data imply that the 
timescale for efficient core-envelope angular momentum transport is substantially less than 
the age of the Sun, but comparable to the lifetimes of the relevant phases (of order 20- 
100 Myr); i.e. solid-body rotation enforced at all times is a poor fit (Keppens et al. 1995; 
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Krishnamurthi et al. 1997; Allain et al. 1998; Sills et al. 2000). Eggenberger et al. (2005) 
found that their diffusion coefficients predicted essentially solid-body rotation at all times, 
which would not be consistent with young open cluster data. We defer a detailed discussion 
of the spindown properties of the revised prescription to a paper in preparation, but we will 
present evidence here that the predicted timescale for core-envelope coupling is longer in the 
new prescription than in the older one. Adopting the revised coefficients therefore improves 
the consistency between theoretical models and the data, by permitting the observationally 
required deviations from solid-body rotation. Combining these two sets of constraints, we 
therefore conclude that some other mechanism is required to explain the slow rotation of the 
solar core, either another magnetic instability or wave-driven momentum transport; however, 
inclusion of the Tayler-Spruit instability may be a promising explanation for the timescale 
of core-envelope coupling. 

As a test, we compare our revised transport coefficients for the Tayler-Spruit dynamo 
to previous results for a solar model. We employ the stellar evolution code described by 
Denissenkov et al. (2006); gravitational settling was not included. We used the Grevesse 
& Noels (1993) mixture of heavy elements and calibrated our model to reproduce the solar 
luminosity (3.85 x lO^^ergs"^) and radius (6.96 x 10^° cm) at the solar age of 4.57 Gyr; this 
procedure yielded Y = 0.273, Z = 0.018, and a mixing length a of 1.75. 

For angular momentum evolution, it is important to specify the initial conditions and 
boundary conditions. Our computations start on the pre-MS deuterium birthline of Palla & 
Stabler (1991) with an initial rotation period of 8 days, close to the mean value for classical 
T Tauri stars. Solid-body rotation was enforced in convective regions at all times. No 
interaction between the protostar and disk was included, implying that the model considered 
corresponds to a rapidly rotating young open cluster star. We adopt the angular momentum 
loss prescription of Krishnamurthi et al. (1997): 



where R and Qg are the star's surface radius and angular velocity, f^crit is the critical velocity 
at which angular momentum loss from the magnetized stellar wind gets saturated. Following 
Krishnamurthi et al. (1997), we adopt ficrit = 10 ^q, where ^q ^ 2.87 x lO^^rads^"*^. The 
parameter K„ is calibrated requiring that our evolved models have Qg = ^q at the solar age 
(-ft'w = 1-1 X lO'^^cm^gs^^ for our unmixed model). 

With these basic ingredients we consider three angular momentum transport cases: a 
model without transport in the radiative interior; a full evolutionary model using the Maeder 
& Meynet (2004) coefficients; and solutions with the revised coefficients employing static 
solar mass and abundance models at different ages. 
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For comparison with the rotation profile obtained by Eggenberger et al. (2005), we used 
the old relations (14- 16), together with the equation of the transport of angular momentum 
by magnetic viscosity: 

In our revised diffusion coefficient estimates the quantities r^c and v^ depend strongly on 
both f2 and g, and there is a steep reduction in the efficiency of angular momentum transport 
in the presence of /i gradients. These features make it numerically more challenging to 
calculate them and q with a reasonable accuracy. However, we can estimate a lower bound 
to the solar rotation curve reliably under the following set of assumptions. First, we assume 
that the timescale for angular momentum transport is less than the solar age in the absence 
of /i gradients. However, there is a threshold value of the /i-gradient above which the new 
coefficients vanish. Assuming that r^c = in (13), the threshold /i-gradient can be estimated 
from the following equation: 

(91nyU 



K = 9V 



dr 



n'q', (33) 



where g is the local gravity. Although formally r]e cannot be less than the magnetic diffusivity 
f], the ratio 1]/K is so small that equation (33) gives a very accurate estimate of the threshold 
/i-gradient. We integrated this condition from the surface to the point where the local 
rotation would be higher than the no-transport case for three different snapshots of the solar 
composition profile (at 150 Myr, 600 Myr, and 4.57 Gyr). 

Our results are summarized in Fig. 3. In the bottom panel, we compare the normalized 
solar rotation as a function of radius obtained from the helioseismic data of Couvidat et 
al. (2003) (dots with errorbars) to a model without angular momentum transport (the long- 
dashed line), a full evolutionary model using the Maeder & Meynet (2004) prescription (solid 
line), and rotation curves derived from static models with /x profiles appropriate for solar 
models at ages of 150 Myr, 600 Myr, and 4.57 Gyr (dotted, short-dashed and dot-dashed 
hues, respectively). In the top panel, we compare the magnetic viscosities that we would 
have obtained for the no-transport rotation curve using our prescription (the solid line) and 
the Maeder & Meynet (2004) prescription (the dashed line). The results in the top panel 
indicate that magnetic angular momentum transport would be highly effective for stellar 
envelopes for both the new and old formulations of the Tayler-Spruit instability. The new 
diffusion coefficients, however, are 2-3 orders of magnitude lower than the previous values 
outside the solar core, and vanish completely in the deep solar interior. 

When compared with the solar data, the marginal stability cases for magnetic transport 
with our diffusion coefficients are clearly well above that permitted by helioseismology. This 
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effect is significant even for composition gradients appropriate for tlie early MS (tlie 150 Myr 
and 600 Myr curves). Tliis mechanism is tlierefore not responsible for the slow rotation of the 
solar core. However, the reduction in the inferred envelope diffusion coefficients may bring 
them into a range compatible with the spindown of young cluster stars. We conclude that, at 
least for the evolution of low mass stars, inclusion of this magnetic transport mechanism may 
have an interesting impact on some observationally important features of models including 
rotation, but that some other effect is still missing in our understanding of the internal solar 
rotation. 

We are grateful to Grant Newsham, Don Terndrup and two anonymous referees for their 
comments and suggestions that have served to improve this paper. We also thank people 
who commented on the first version of our paper in astro-ph/0604045. Tony Piro brought 
our attention to subtle fiaws in our original argument based on the dispersion relation. 
Andre Maeder mentioned the relationship between the energy of mixing and the energy of 
differential rotation, and Henk Spruit provided helpful insight into the dispersion relation. 
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Fig. 1. — Length scales of the displacement caused by the Tayler-Spruit instability. 
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Fig. 2. — An illustration of the motion of a magnetic ring of the radius r = 2 (in arbitrary 
units) constrained by the Coriolis force (eqs. 27). The large circle is the equator of the 
sphere of the radius r = 10 inside the star. The coordinate plane XOY is perpendicular to 
the rotation axis. The smaller circles show consecutive positions of the ring at the phases 
2Qt = 0, tt/2, it, and 37r/2. The ring starts its motion from the pole with the velocity 
V = {va, 0}, where Va = ruA. The plotted case corresponds to the ratio Ua/Q = 0.1. 
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Fig. 3. — Panel a: the old {long- dashed curve) and new {solid curve) magnetic viscosities 
for the no-transport rotation profile {long-dashed curve in panel b). Panel b: solid curve — 
rotation profile from our full evolutionary computations obtained as a solution of eq. (32) 
using the Maeder & Meynet (2004) prescription for Ug] dotted, short- dashed and dot- dashed 
curves — rotation profiles for our revised u^ estimated with equation (33) using A^^ from 
models of ages 1.5 x 10^, 6 x 10^, and 4.57 x 10^ yr, respectively; dots are helioseismic data 
from Couvidat et al. (2003). 



